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Abstract 

The running of renormalized quark masses is computed in lattice QCD 
with two flavors of massless 0(a) improved Wilson quarks. The regulariza- 
tion and flavor independent factor that relates running quark masses to the 
renormalization group invariant ones is evaluated in the Schrodinger Func- 
tional scheme. Using existing data for the scale tq and the pseudoscalar 
meson masses, we define a reference quark mass in QCD with two degen- 
erate quark flavors. We then compute the renormalization group invariant 
reference quark mass at three different lattice spacings. Our estimate for 
the continuum value is converted to the strange quark mass with the help 
of chiral perturbation theory. 



1 



1 Introduction 



Lattice QCD provides a definition of quark masses from first principles. For 
example given a hadronic input accessible to experiments, such as the K-meson 
mass TTiK and decay constant Fk, it is possible to compute the strange quark mass 
on the lattice. A convenient quantity to consider is the renormalization group 
invariant (RGI) strange quark mass, which is independent of the renormalization 
scheme if the renormalization conditions are imposed at zero quark mass [1] . 

If one looks at the most recent lattice results for the strange quark mass, from 
simulations with two [2-4] and two plus one [5,6] dynamical quark- flavours, there 
is a large spread of values ranging from 68MeV to 132 MeV. There are several 
sources of systematic errors in these computations, such as the use of perturbative 
renormalization (except for [4] , where non-perturbative renormalization is done in 
the RI-MOM scheme) and the values of the lattice spacing affordable nowadays. 

To perform a completely controlled computation we start from the bare 
PC AC mass mi{go) for a given quark flavor i and determine the RGI mass 

Mi = Zuigo) mi{go) . (1.1) 

The bare coupling qq is in a one-to-one relation to the lattice spacing and the 
continuum limit of Mj exists (and should be taken). As an intermediate step we 
first define the running mass 

MiA - ^r^,m{go) . (1.2) 

where rn is non-perturbatively defined in the Schrodinger Functional renormal- 
ization scheme and hence is well-defined also for low energies This corresponds 
to splitting Zm into two factors 

which will be computed with full non-perturbative precision following the strategy 
of [7]. Note that all renormalization factors are flavor-independent and we hence 
omit the subscript i in mass ratios. In the flrst factor of the above splitting the 
universal continuum limit is understood to have been taken. Its computation is 
the main objective of this paper. This result can then be used for any action 
and only the second factor needs to be redetermined. This becomes tractable by 
making the universal factor available for /i in the range of hadronic energies of 
0(1 — 2 GeV). As for the running couphng [8] our method thus avoids the need 
to treat a multi-scale problem in a large volume in this step. 
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First results on the /i-dependence of m in the N{ = 2 theory have already 
appeared in Refs. [9-11]. Beyond a finalization of these results, we here compute 
the second factor — and hence Zm — using non-perturbative improvement [12] 
and renormalization [13] of the axial current in the theory with two flavors of 0(a) 
improved Wilson quarks and plaquette gauge action. 

As an application, one wants to determine the light (up, down, strange) 
quark masses.^ For practical reasons this is at present not yet possible for us 
in the most straight-forward way by simulations at the physical parameters. We 
need a hierarchy of additional approximations to compute for instance the strange 
quark mass. 

First of all our simulation algorithm is at the moment still restricted to pairs 
of degenerate flavours and we include one such pair (A^f = 2). In the large volume 
simulations needed to determine mi{go) in Eq. (1.1) our masses can at present 
not be taken to values small enough for the up- and down-flavors, for instance 
by tuning the pseudoscalar states to the physical pion mass. Instead we shall 
determine the RGI quark mass Mref that is associated with a "Kaon" made from 
two degenerate flavours. 

The analogous computation of Mref has been performed in the quenched 
A^f = theory [17] with a result that agrees within errors with the new one. We 
therefore assume that at the present level of accuracy it also applies to a hypo- 
thetical QCD with three degenerate flavors. A relation between this model and 
real QCD is flnally established by chiral perturbation theory [18,19] supplemented 
with some knowledge [19-21] of the phenomenologically inaccessible [22] low en- 
ergy constants. The conclusion of [19] is that the ratio of quark masses is close to 
the one given at lowest order in chiral perturbation theory. Given this "fact" — 
but keeping in mind that it should be scrutinized in future lattice QCD compu- 
tations — it is then sufficient to compute any one quark mass from lattice QCD. 
In particular we can connect our Mref with the strange mass Mg. The result that 
corresponds to the Gell-Mann-Oakes-Renner formula [23] is Mg — 48/25Mref. 

Finally, the conversion of the RGI mass to the conventionally cited MS mass 
at 2 GeV renormalization scale is of course based on perturbation theory, which 
does however look very well behaved, see Table 6. 

^Later, when lattice spacings are small enough, the charm quark mass can be determined 
[14, 15]. Also the beauty quark mass computed in HQET through the strategy of [16], is based 
on Eq. (1.1) in QCD. 
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2 The renormalization scheme 



QCD is a theory which has as free parameters the bare gauge couphng go and 
the bare quark masses rrii , i = 1,...,N{. The hadronic scales hke Fk or 
are connected to the perturbative high energy regime of QCD via the running of 
renormahzed couphngs 

fill) = ZgQ^, mi{n) = Zmrui, (2.1) 

where // is the renormahzation scale. In the following we assume that the renor- 
malization conditions are imposed at zero quark mass (mass-independent schemes) 
[1]. Introducing a regularization prescription, e.g. a lattice spacing a, the renor- 
malization factors are functions of go and a/i 

^9 = Zg{go, aix) , Zm = Z^{go, a/i) . (2.2) 

In renormahzed quantities the regulator can be removed, e.g. the continuum limit 
a — > can be taken, yielding a finite result. The advantage of mass-independent 
renormalization schemes is that in all such schemes the ratios of renormahzed 
quark masses for different flavors are scale and scheme independent constants 
[19,24]. 

The running of the renormahzed couplings is described by the renormalization 
group equations (RGE) 

^% ^ ' ^ 'T{9)rni . (2.3) 

The P and r functions are non-perturbatively defined if this is true for g and mj. 
Their perturbative expansions are 

m ~ -f{bo + bif + b2f + ...}, (2.4) 

r{g) ^ -f{do + drf + ...]. (2.5) 



The coefficients 



are scheme independent. A physical quantity P is a quantity for which the total 
dependence on the renormalization scale vanishes, i.e. it is a renormalization 
group invariant (RGI) 

/i^P(/x,^,{m,}) = 0. (2.7) 
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Examples are the A-parameter and the RGI quark masses 
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m 

dp 
box 



hlx 



(2.9) 



where g = g{fi) and rfii = rni{fi). The A-parameter and the RGI quark masses 
are defined independent of perturbation theory and their connections between 
different mass independent renormahzation schemes can be given in a simple and 
exact way [24]. In particular the RGI quark masses Mj are scheme independent. 
Any physical quantity P can be considered to be a function of A and Mj, i.e. 
there exist a function P such that [24] 

P{l,,g,{mi}) = P{A,{Mi}). (2.10) 

For this reason A and Mj seem preferable as the fundamental parameters of QCD. 



2.1 Quark masses in the Schrodinger Functional 

In QCD a renormalized mass is defined through the partially conserved axial 
current (PCAC) relation, which involves the renormalized axial current (Ar)^(x) 
and the renormalized pseudoscalar density Pr{x) 

5m(^r)m = {m + mj)P^, (2.11) 
iA^)^ix) = ZA^,(a;)7^75^,(x), (2.12) 
Pn{x) = Zp^,(x)75^,(a;). (2.13) 

The renormahzation constant Za can be calculated non-perturbatively, using chi- 
ral Ward identities and does not depend on the renormahzation scale. 

The renormahzation constant Zp can be conveniently determined in the 
Schrodinger Functional (SF) renormahzation scheme [25,26]. There QCD is for- 
mulated in a finite box of spatial size L and temporal extent T. The fields are sub- 
ject to Dirichlet boundary conditions in time, which provide an infrared cutoff to 
the frequency spectrum of quarks and gluons. This allows to perform simulations 
at zero quark mass and thus to use the SF as a mass-independent renormahzation 
scheme. Our renormahzation scheme is further specified by setting T = L (see 
below). The renormahzation conditions are then naturally imposed at the scale 
= 1/L. 

The presence of boundary values for the fields in the SF formulation of a field 
theory requires in general additional (compared to the case without boundaries) 
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xo = 



■ ^ 

LxLxL 



Figure 1: The Schrodinger Functional. 

renormalization [27,28]. The renormahzabihty of QCD with SF boundaries has 
been studied in Refs. [29-32] where it was shown that no additional counterterms 
are needed except for one boundary term which amounts to a rescaling of the 
boundary values of the fermion fields by a logarithmically divergent factor. 

The renormalization condition for Zp that we employ is discussed in Refs. 
[7,24,25]. It uses a correlation function /p(a;o), which is a matrix element of the 
pseudoscalar density inserted at time distance xq from a pseudoscalar boundary 
state, the other boundary state having vacuum quantum numbers. To cancel 
the multiplicative renormalization of the boundary quark fields the boundary-to- 
boundary correlation function /i is used. The renormalization constant Zp is then 
defined through 

^P = 7T7^ m, = 0, z = l,...,iVf. (2.14) 

The correlation functions are schematically represented in Fig. 1. They are com- 
puted at zero quark masses mj. The definition Eq. (2.14) is such that Zp = 1 at 
tree level of perturbation theory. The renormalization condition for Zp is further 
specified by setting 

T = L, C = C' = 0, e = 0.5, (2.15) 

where C , C are the boundary gauge fields in the Lie algebra at xq = and xq = T 
and 6 the parameter controlling the spatial boundary conditions of the fermion 
fields. For more details of the calculation we refer to [7,24]. 

A rigorous definition of the renormalized mass can be given in the lattice 
regularization of QCD. In our case we work with Nf = 2 mass-degenerate flavors 
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of 0(a) improved [33,34] Wilson fermions. The massless theory is defined in the 
bare parameter space along the line k = Hcido) where the PC AC mass 



, , l{d^ + do)fA{xo) + CAad^dofp{xo) 



2/p(a:o) 



(2.16) 

XQ=T/2 



vanishes. Here, Oq and Oq are the forward and backward lattice derivatives, re- 
spectively. The 0(a) improvement coefficient Ca of the axial current has been 
computed non-perturbatively in Ref. [12]. The correlation function f\{xo) of the 
axial current Aq is defined analogously to fp{xo). The renormahzed PC AC mass 
at the scale /i — 1/L can then be written as 

m{ii) = lim Z^{go,aii)m{go,K) , Z^{go,aiJ,) = ^a(5'o) ^ ^2.17) 

o-»0 u=g^{L) z-p[go,L/a) 

where the renormahzed gauge coupling g^{L) [8] is kept fixed. Since the coupling 
runs with the box size L, keeping g^ fixed means keeping the renormalization scale 
/J, fixed. The renormalization factor of the axial current has been determined 
non-perturbatively in Ref. [13]. The precise definition of Zp (differing at order 
fi-om Eq. (2.14)) is 

Zp{go,L/a) = c / at i^ = Kc, (2.18) 

where the factor c{L/a) is chosen such that Zp{0,L/a) — 1 and can be found 
in [7]. 

In order that the continuum limit is reached with cutoff effects strictly pro- 
portional to a^, the 0(a) improvement factor {l + {bA — bp)amq) should be included 
in Eq. (2.17) [33]. Here am^ = (l//t— 1/kc)/2 is the bare subtracted quark mass. 
Results in perturbation theory [35] and in the quenched approximation [36] show 
that the difference of improvement coefficients 6a — bp is small. The quark masses 
in our simulations will also be relatively small and we expect corrections due to 
— bp at the per mille level, which we neglect. 



3 The running of the mass in the SF— scheme 



The running of the renormahzed quark mass m(//.) in the SF scheme as specified 
in Section 2.1 with Nf — 2 mass-degenerate fiavors can be computed on the lattice 
from the step scaling function of the renormalization factor Zp extrapolated to 
the continuum 

Zp{go: 2L/a) 



(Jp{u) = lim Sp(m, a/L) , 'Ep(u,a/L) 



Zp(go, L/a) 



u=g^L) 



(3.1) 
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u 



ap{u) 



0.9793 
1.1814 
1.5031 
2.0142 
2.4792 
3.3340 



0.9654(9)(11) 
0.9527(11)(6) 

0.9413(16)(2) 
0.9174(16)(24) 
0.8871(23)(18) 
0.8384(35)(12) 



2.16 
0.47 
0.01 

3.58 
0.54 
0.20 



Table 1: Continuum extrapolations of Ep fitting the L/a — 8 and L/a — 12 data 
to a constant. The first error is statisticaL The second error is the difference 
between the fit and the L/a — 8 results and will be added linearly as a systematic 
error. 

Prom Eq. (2.17) it follows immediately that 



i.e. the step scaling function ap{u) describes the running of the renormalized quark 
mass. We computed Sp(M,a/L) at six values of the renormalized coupling u cor- 
responding approximately to a range of box sizes of the order L = 10~^ fm . . . 1 fm 
(or equivalently /j, of the order 100 GeV ... 1 GeV). At each value of u we simu- 
lated three lattice resolutions L/a — 6, 8, 12 and the results for Zp and Ep are 
summarized in Table 7 in Appendix A. For the extrapolation to the continuum, 
we fitted to a constant the two values of Sp on the finer lattices, separately for 
each coupling u. We then added linearly the difference between the fit and the 
L/a = 8 result as a systematic error. The continuum estimates can be seen in 
Fig. 2. Our data do not show any significant dependence on the lattice spac- 
ing, as we could verify by trying different extrapolations (quadratic, linear in a). 
This statement is based on the statistical accuracy that we could achieve. We re- 
mark that also in the quenched approximation the cutoff effects were found to be 
small [7] and there Sp was computed at an even finer lattice resolution L/a = 16. 
The continuum values of crp{u) arc summarized in Table 1. In the last column 
we list the divided by the number of degrees of freedom ndf of the fit. Their 
average is close to the expected value of one. 

In perturbation theory the step scaling function ap{u) has an expansion 
ap{u) = 1 — ln(2)doti + 0{u^). In Fig. 3 our non-perturbative data for ap{u) are 
conveniently plotted for comparison with perturbation theory. We relate ap{u) to 



ap{u) 



m{fi) 



for fj, — 1/L , 



(3.2) 



m(/./2) 



8 



0.98 
0.96 
0.94 
^ 0.92 
5 0.9 



0.88 
0.86 
0.84 
0.82 







- - ^ — 


© 




- © 




— © 




- -®- - 


— a> — 




<^-- 


® 




© 



0.01 



0.02 

(a/L)2 



0.03 



0.04 



Figure 2: Continuum extrapolations of Sp. The L/a = 6 data have been excluded 
from the fits. For the third smallest coupling, the two points at L/a = 6 refer to 
1-loop and 2-loop values for the boundary improvement coefficient q [29,37]. 



r and P using Eq. (2.9) 



(TplM) 



u 



a{u) 



do/{2bo) 



exp 



r{x) _ do_ 
(3{x) box 



(3.3) 



where cr{u) is the step scaling function of the coupling and is determined by [8] 



-21n(2) 



a(u) 



dx 



^/x(3{^yx) ' 

Using for the r-function the 2-loop expression with [24] 

di = 4(0.0271 + O.OlOSiVf) 
and for the /3-function the 3-loop expression with (see [8]) 

62 = (0.483 - 0.275iVf + 0.0361A/'f2 - 0.00175iVf^)/(47r) 



(3.4) 



(3.5) 



(3.6) 



we get from Eq. (3.3) and Eq. (3.4) the perturbative curve shown in Fig. 3. Our 
non-perturbative data do not show any significant deviation from the perturbative 
estimates. 



9 




Figure 3: Comparison of the non-perturbative data for the step scahng function 
(Jp{u) with perturbation theory. 2/3-loop refers to the 2-loop r-function and 
3-loop /^-function, analogously 1/2-loop. 



We now take the continuum values for crp{u) in Table 1 (with added statistical 
and systematic errors) and for the step scaling function of the coupling a{u) the 
results from fits to constants in Table 4 of Ref. [8]. We solve the following joint 
recursion to evolve coupling and mass from a low energy scale 1 / Lmax defined by 

uo = ^'(L^ax) = 4.61 (3.7) 

to the higher scales l/Lk, k = 0,1, ... ,8 (with Lq = Lmax) 

I = ^'^^"^^^^ = ^-^^ ^ Uk = f{L,) , L, = 2-'=L_ , (3.8) 
We interpolate the values of (y{u) and ap{u) through a polynomial ansatz 

a{u) = U + SqU^ + SiU^ + S2M^ + S3U^ + SiU^ , (3.10) 

(Tp(n) = 1 + pqu + piu'^ + p2U^ , (3-11) 

where the coefficients sq, si [8] and po = — ln(2)(io are fixed to their perturbative 
values. The coefficients S2, S3, S4 and pi, p2 are here fit parameters. The errors 
of the recursion coefficients are computed by error propagation. 



10 



k Uk M/m(l/L^ax) M/m{l/L^^) 
2/3-loop 1/2-loop 






4.61 


1.274 


1.267 


1 


3.032(16) 


1.296(6) 


1.292 


2 


2.341(21) 


1.295(10) 


1.292 


3 


1.918(20) 


1.294(13) 


1.292 


4 


1.628(16) 


1.294(14) 


1.292 


5 


1.414(14) 


1.295(15) 


1.293 


6 


1.251(12) 


1.297(16) 


1.295 


7 


1.121(10) 


1.298(17) 


1.297 


8 


1.017(10) 


1.299(17) 


1.298 



Table 2: Values for M/m{l/L^^) from Eq. (3.12). 



Using the coefficients Wk in Eq. (3.9) we compute 



- w,'———, (3.12) 



m(l/L^ax) ^ m(l/Lfe) ' 

where the factor M/rn{l/ L^) is calculated from Eq. (2.9) with — by em- 
ploying the perturbative expressions for the r- and /3-functions at 2- respectively 
3-loop order. The results arc shown in Table 2. They have a remarkable stability 
in the coupling^ Uk and we take A; = 6 as our result 

^ -- 1.297(16) at ji = l/L^^. (3.13) 



m{n) 

We emphasize that from Eq. (2.9) it is evident that a-p{u) and M/m{n) are flavor 
independent. Moreover, since the continuum limit has been taken any regulariza- 
tion dependence has been removed from the result Eq. (3.13). 

Finally, in Fig. 4 we plot the non-perturbative running of the renormalized 
mass. For fi/A = l/(LfcA), = 0, 1, . . . , 8, we plot the points rn{l/ Lk) / M ob- 
tained from Eq. (3.12) using the result Eq. (3.13). The physical scale A is here 
implicitly determined through In(ALmax) = —1.298(58) obtained from the recur- 
sion Eq. (3.8). In the plot we neglect the overall uncertainties of m{l/ Lmia)/M 
and In(ALinax), since they would simply change the scales on the plot axes. The 

^The deviation in the case A; = is due to the difference between the perturbative and the 
non-perturbative values of a{u) at large u (see [8]). 
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1 10 100 1000 

Figure 4: The non-perturbative running of m. 2/3-loop refers to the 2-loop 
r-function and 3-loop /3-function, analogously 1/2-loop. 

errors of the points in Fig. 4 come from the coefficients Wk- Together with the 
non-perturbative points we show the perturbative curves that are obtained from 
Eq. (2.8) and Eq. (2.9) by using the perturbative expressions for the r- and 
/^-functions at 1-, 2- respectively 2-, 3-loop order. The non-perturbative and 
perturbative running are very close down to the smallest energies that were ac- 
cessible in our simulations. We remark that this statement explicitly refers to the 
special SF renormalization scheme considered here. 

4 Estimate of the strange quark mass 

4.1 Complete Zm for the improved Wilson discretization 

We now derive the second factor in Eq. (1.3) for a few values of the lattice spacing 
or respectively the bare coupling. As emphasized before, this contribution is 
non-universal and in the form given it will be valid only for our action of non- 
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13 


K 


L/a 






5.20 


0.13600 


4 


3.65(3) 


0.55188(51) 


5.20 


0.13600 


6 


4.61(4) 


0.47876(47) 


5.29 


0.13641 


4 


3.394(17) 


0.57643(50) 


5.29 


0.13641 


6 


4.279(37) 


0.50739(58) 


5.29 


0.13641 


8 


5.65(9) 


0.45799(75) 


5.40 


0.13669 


4 


3.188(24) 


0.60269(65) 


5.40 


0.13669 


6 


3.861(34) 


0.53681(54) 


5.40 


0.13669 


8 


4.747(63) 


0.49116(66) 



Table 3: Results for Zp, q set to 2-loop value. The values of are from [8]. The 
hopping parameters k are set to the critical ones (kc) of [4]. 

perturbativcly improved Wilson fermions with plaquette gauge action and Cgw 
as specified in [34]. For and ca we employ the values recently given and 
parameterized in [13] and [12]. It remains to compute Zp((7o, -^^max/a) for the 
desired values of the bare coupling g^^ here given by /5 = 5.2, 5.29, 5.4. The scale 
Lmax is fixed by ^^(Lmax) = 4.61, where the universal factor of Zm is now known, 
Eq. (3.13). Our basis are the simulation results summarized in Table 3. 

While the simulation at the largest bare coupling is exactly at the target value 
for (?^, the two other series of simulations require a slight interpolation. This has 
been done using a fit ansatz motivated by Eq. (2.9) 

In(Zp) = Ci + C2ln(^') (4.1) 

to interpolate Zp between two values of cp' straddling 4.61. The fit takes into 
account the (independent) errors of both Zp and cp' . The fit error is then aug- 
mented by the difference between the fit result from Eq. (4.1) and the result from 
a simple two point linear interpolation in g^ . The values of the coefficient C2 in the 
fit Eq. (4.1) are found to be -0.369(25) at /3 = 5.29 and -0.430(34) at /3 = 5.40, 
which are not far from — (io/(26o) = —0.4138. 

The resulting numbers for Zp and Zm are summarized in Table 4. The first 
error for Zm comes from the error of the factor Zp^jZ^. The second error is the 
1.2% uncertainty in the universal factor M/m and should be added in quadrature 
to the quark mass error after the continuum limit has eventually been taken. 
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(3 Zp Zyi 



5.20 0.47876(47) 1.935(33) (24) 
5.29 0.4936(34) 1.979(25) (24) 
5.40 0.4974(33) 2.001(29) (25) 

Table 4: Results for Zp and finally Zm for three coupling values. The action to 
which it refers is detailed in the text. 

4.2 The reference quark mass 

As announced in the introduction we next compute the reference quark mass M^f 
producing a pseudoscalar with the mass of the Kaon in our simulated two-flavour 
theory. The hopping parameter K^ei is tuned to keep the pseudoscalar mass mps 
in the relation 

(ro(«;c)"^Ps(«ref))^ = (romx)^ = 1.5736, (4.2) 

where rg is the scale extracted from the static quark potential [38]. This value 
corresponds to the K-meson mass = (495 MeV)^ for ro = 0.5 fm [17]. Using 
the data for ro(«;), mps(fi:) and available in Ref. [4] and the extrapolations to 
the chiral limit ro(Kc) of Ref- [8]; we can safely perform a slight extrapolation to 
the value k — K^ei defined in Eq. (4.2). We employ the following fit ansatz for 
romps as a function of the bare subtracted quark mass amq 

(ro(fi;c)mps(K))2 

= ei + 62 • arria , (4.3) 

with fit coefficients Ci and 62- The results of the ffis are /tfcf = 0.135680(30) at 
p = 5.20, Kref = 0.136018(27) at /3 = 5.29 and n^ei = 0.136293(25) at /? = 5.40. 

The PCAC masses OTTT^ef have been computed in simulations at — Href 
and the results are summarized in Table 5. In these simulations the parameter 
6 in Eq. (2.15) has been set to zero. Barring cutoff effects the PCAC mass is 
independent of the time Xq at which the right hand side of Eq. (2.16) is evaluated. 
Therefore in Table 5 we average over a range ti : t2 of xq values around xq = 
L/2, where the PCAC mass has a plateau. We keep ti — t2 roughly constant 
in physical units. We see a significant gain in statistical precision due to the 
averaging compared to taking only the time Xq — L/2. The error analysis of 
derived observables like the (averaged) PCAC mass and Zp has been done with 
the method of Ref. [39]. 
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p 




L/a 


ti/a : t2/a 


amref 


Mref[MeV] 


M,ef[McV] 














using Z^"" 


5.20 


0.135680 


16 


7:9 


0.01410(30) 


58.7(3.2) 


69.5(3.7) 


5.29 


0.136018 


16 


6:10 


0.01352(28) 


63.5(3.2) 


69.1(3.6) 


5.40 


0.136293 


24 


10:14 


0.01300(18) 


72.0(2.7) 


76.2(2.7) 



Table 5: Results for the PCAC masses and the RGI reference quark masses, the 
latter converted to MeV assuming tq — 0.5 fm. Two definitions of Za are used [13] . 



The cutoff effects in mref also depend on the volume. In fact, at a lattice 
spacing of 0.1 fm this dependence is rather strong [40]. Thus one has to define 
the PCAC mass for a fixed volume in order to ensure that cutoff effects disappear 
smoothly as O(a^). Furthermore it is preferable to choose L relatively large. We 
choose L > 1.5 fm, where the volume dependence can be neglected. Indeed, we 
compared results from a simulation at = 5.2, k, = 0.1355, Csw = 2.02 on a 
L/a = 16 lattice with those from the JLQCD collaboration"^ obtained at the same 
parameters on a 20"^ x 48 lattice [3]. We could check that at the level of 1.6% 
statistical precision in our simulation, the PCAC masses agree. We conclude that 
our volume is large enough at /3 = 5.2 for the volume dependence to be negligible. 
Therefore at the higher (3 values we choose approximately matched (or larger) 
physical volumes. 

The second to last column of Table 5 shows our results for the reference RGI 
quark mass 

M,ef = ZMrUrei. (4.4) 

We do see large cutoff effects. At /? = 5.2 similarly large cutoff effects have been 
observed in other quantities as well [13,40,41]. Moreover the lattice spacing in 
our simulations changes by 30% only when going from P — 5.2 to P — 5.4. We 

therefore do not attempt (and discourage) any elaborate continuum extrapolation 
based on this range of lattice spacings. Instead we use a conservative estimate of 
the continuum value of Mref by taking its value at our largest (3 = 5.4 with the 
difference to the value at /9 = 5.2 added as systematic error 

Mref = 72(3) (13) MeV. (4.5) 

This result makes it clear that at present the systematic error is dominating over 
the statistical one. 

^We are grateful to Takashi Kaneko for providing us with data for this comparison. 
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The last column of Table 5 shows the results for an alternative estimate of 
Mref, which differs only by the use of the renormalization factor Z^°^ [13]. We 
observe that the difference in Mref between the two choices for is compatible 
with an a^-behavior. Actually, the data obtained with have a weaker a- 
dependence but this effect is essentially due to the change at the coarsest lattice 
spacing {(3 = 5.2). 

4.3 The strange quark mass 

To make contact with physics we now assume that Eq. (4.5) also holds within 
errors in a theory with three degenerate flavors, which can then be related to 
the strange quark in QCD by chiral perturbation theory This assumption is 
supported by the fact that the value of M^f Eq. (4.5) is the same within errors 
as in the A^f = theory [17]. 

At lowest order in chiral perturbation theory, disregarding the electromag- 
netic interaction, the formula [18,23,42] 

= ^(^K++^Ko) = (M + M,)5rgi, (4.6) 

where M = 1/2 (Mu + Ma), holds. Here M^, Md, Ms are the up, down, strange 
RGI quark masses and -Brgi is a constant of the chiral Lagrangian. Eq. (4.6) 
imphes for degenerate quarks of RGI mass Mref (defined according to Eq. (4.2)) 

m| = 2Mref5RGI. (4.7) 

Therefore the relation Mref — {M + Ms)/2 holds at lowest order in chiral pertur- 
bation theory and using Mg/M = 24.4(1.5) [19] gives 

Ms ^ 48/25 iV/ref. (4.8) 

Corrections to Eq. (4.8) are expected to be small in chiral perturbation theory 
[17,20,43,44] and anyway below the accuracy that we will reach here for our result. 
In addition in the quenched approximation it was found that the dependence of 
pseudoscalar masses rrips, at fixed average quark mass, on the difference of quark 
masses is rather small [17]. We therefore assume Eq. (4.8) to hold also with 
dynamical quarks. 

At this point we are ready to give an estimate of the continuum value of the 
RGI strange quark mass 

Ms = 138(5) (26) MeV, (4.9) 
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//[GeV] m^^{fx)/M 





2-loop 


3-loop 


4-loop 


1.0 


0.7979 


0.8364 


0.8469 


2.0 


0.6824 


0.6984 


0.7013 


4.0 


0.6110 


0.6197 


0.6209 


8.0 


0.5608 


0.5662 


0.5668 


90.0 


0.4571 


0.4588 


0.4589 



Table 6: Factors to convert the renormalization group invariant mass into the MS 
scheme at scale /x, for A^^ = 245(32) [8] and A^f = 2. 

by combining Eq. (4.8) with Eq. (4.5). Equivalently to the determination of M^ef, 
from Eq. (4.7) and Eq. (4.5) we get 

Bji^Gi = 1.70(38) GeV. (4.10) 

The renormahzed strange quark mass in the MS scheme at the renormahza- 
tion scale jj, is obtained by multiplying Mg with the conversion factor m'^S(/i)/M. 
The latter is computed perturbatively by numerical integration of Eq. (2.8) and 
Eq. (2.9) and is listed in Table 6 for different choices of /i. The n = 2, 3, 4-loop 
approximations of the f3 and r functions in the MS scheme (for Nf = 2) are used. 
The coefficients at 4-loop have been computed in [45-48]. To set n in physical 
units the result = 245(32) of Ref. [8] is used. The uncertainty in the A pa- 
rameter translates into a 2.8% uncertainty in m'^^(//)/M at // = 2 GeV and 4.9% 
at n = 1 GeV (with 4-loop evolution) . Taking this into account our estimate for 
the MS strange quark mass is 

mf^(/x) = 97(22) MeV at /x = 2GeV. (4.11) 
5 Conclusions and outlook 

We have presented a fully non-perturbative renormalization of the quark mass 
in two flavor QCD, with renormalization conditions specified at zero quark mass 
in the Schrodinger Functional scheme. Simulations were performed with 0(a) 
improved Wilson quarks. Our main results are the running of the quark mass in 
Fig. 4 and the factor M/m(//) relating the quark mass at a specified low energy 
scale II with the RGI quark mass in the continuum limit, Eq. (3.13). 
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Figure 5: Summary of the strange quark mass data from lattice simulations. In 
the legend the discretizations used are indicated in the form gauge action/fermion 
action. The dictionary reads: W: Wilson gauge action; I: Iwasaki gauge action; 
LW: 1-loop tadpole improved Liischer-Weisz gauge action; CW: Wilson-clover 
fermion action; KS: Asqtad staggered fermion action. The dotted lines represent 
the quenched result [17]. 

In order to obtain the renormalized strange quark mass in physical units an 
appropriate hadronic scheme has to be defined. By using existing data on the 
pseudoscalar masses and the scale vq from Ref. [4] we were able to determine the 
RGI reference quark mass Mref defined through Eq. (4.4) and Eq. (4.2), at three 
lattice spacings in the approximate range 0.092 . . . 0.071 fm. Even in this small 
range we see distinct cutoff effects. At the largest lattice spacing large cutoff 
effects have been observed elsewhere [13,40,41]. These facts make it impossible 
to perform a systematic continuum extrapolation. Nevertheless we can give a 
conservative estimate of the continuum value in Eq. (4.5). We emphasize that the 
error in Eq. (4.5) is dominated by the systematic uncertainty of this step. 

Using a number of additional reasonable assumptions we convert Mref to the 
strange quark mass in the MS scheme Eq. (4.11). A definite clarification that all 
of these lead to negligible errors will require some future work. 

In Fig. 5 we summarize the most recent results for the strange quark mass 
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from lattice simulations with two [2-4] and two plus one [5, 6] dynamical quarks. 
The physical K-meson mass is always used as input although in some cases this 
is employed in a partially quenched setup where valence and sea quark masses 
differ, while we compute Mj-cf and then use chiral perturbation theory to connect 
to the physical theory. In the legend we emphasize whether nonperturbative or 
perturbative renormalization has been used as well as the adopted discretizations. 
It appears that perturbative renormalization leads to rather small values for the 
strange quark mass.^ We plot the results of our present work using two renormal- 
ization conditions for Za that differ at O(a^). The red triangles are obtained with 
[13] and are slightly displaced for clarity. This comparison gives a flavor of 
the quite large cutoff effects at /3 = 5.2, our coarsest lattice spacing. The dotted 
lines in Fig. 5 mark the quenched result of Ref. [17]. Given the present status, 
illustrated in Fig. 5, it appears hard to claim a definite dependence of mg on the 
number of dynamical fermions, even between A'^f = and = 2. 

Our present result for the renormalized strange quark mass should be im- 
proved by simulating at a finer lattice spacing. We are confident that this will 
be possible in the near future. There are promising improvements in algorithms 
[49-51], and a new machine apeNEXT is becoming available to us [52,53]. More- 
over the determination of the lattice spacing can be improved by, for example, 
computing the renormalized strange quark mass in units of the K-meson decay 
constant Fk. There are first indications that the SF is an efficient setup for this 
low-energy computation (see also [54]). 
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Table 7: Results for the step scaling function Ep. 



A Simulation results for Zp 

In Table 7 we collect the bare parameters and results of our simulations to compute 
Zp. Simulations on. L/a and 2L/a lattices are required to extract the step scaling 
function Ep Eq. (3.1). At the three lowest couplings g'^{L) simulations have been 
performed using the 1-loop value of Ct [29], except for L/a = 6, f3 — 7.5420 and 
L/a — 8, P — 7.7206. For the latter parameters and the larger couphngs the 

^Note that the use of lowest order chiral perturbation theory in our work is not hkely to 
be a significant source of difference to the other computations, since they do not report large 
deviations from lowest order chiral perturbation theory, see for example Ref . [3] . 
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Figure 6: Reversibility violations in the HMC as a function of the trajectory length 
r used in the reversibility tests. Empty symbols represent the average and filled 
symbols the maximal value of the Hamiltonian violation | Aif (cycle) | [11]. 



2-loop value of Ct [37] has been used. At the third lowest coupling u ^ 1.5 we 
checked at L/a = 6 that there is no significant difference in Sp using the 1- or 
2-loop value for Ct, as it is shown in Fig. 2. 

To the statistical error of Sp we added in quadrature the error due to the 
uncertainty in the coupling u. The latter was estimated using the 1-loop result 
— ln(2)(io for the derivative of Sp with respect to u. This correction is tiny, it 
increases the errors of Sp by at most 5% at the largest coupling and lattice. 



B About the algorithm 



Our simulations have been performed using the Hybrid Monte Carlo algorithm 
with two pseudo-fermion fields as proposed by M. Hasenbusch [11,55,56]. Here 
we extend the study on algorithmic precision presented in Ref. [11]. 

For our largest lattice volume L/a = 24 at /5 = 5.4 we investigated the issue of 
reversibility. We looked at the quantity | Aif (cycle) | which measures the difference 
of the Hamiltonian for a cyclic trajectory [11]. In Fig. 6 we present data for the 
average | Aif (cycle) | (empty symbols) and the maximum of |Aif(cycle)| (filled 
symbols) on three sets of configurations. Sticking to single precision arithmetic 
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we found for a trajectory length r = 0.5 (that we use in the production run) an 
average value for |Ai7 (cycle) | (black triangles) that is larger by almost a factor five 
with respect to the simulations at /3 = 6.118 reported in Ref. [11] (blue circles). 
We therefore repeated the simulation at /3 = 5.4 in double precision arithmetic. At 
the same time we reduced the parameter e^, defined® as the requested accuracy in 
the conjugate gradient iteration, from = 10~^^ (as used with single precision) 
to = 10~^^. The results for double precision are the red squares in Fig. 6. 
By comparing single with double precision, both the average and the maximal 
value of |Ai?(cycle)| decrease by two orders of magnitude for double precision 
and these values are well below the values observed in Ref. [11]. Our statistics 
for the reversibility checks consists of 18 configurations analyzed at /9 = 5.4 (both 
for single and double precision) and of 30 configurations analyzed at j3 = 6.118. 
On a subset of the (3 = 5.4 configurations we checked that within single precision 
by reducing to 10~^^ (which one would expect to be the lower bound in single 
precision) the roundoff error in the Hamiltonian computation and consequently 
the value of | Aif (cycle) | remain unchanged. 

Finally, we emphasize that we found no significant difference in the observ- 
ables we computed at /9 = 5.4 between the single and double precision simulations. 
The number quoted in Table 5 for mref is from the simulation in double precision. 
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